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Needs ["Laurent Functions * " ] 



RuleDelayed: :rhs : Pattern t_ appears on the right-hand side of rule 

PhaseAngle[L_] [t_] { PhaseAngle [L] [t_] =Module[{xl, x2 , x3 , x4 , x5 , x6}. 

Needs [ " LaurentNotationTes t ^ " ] 

Needs: rnocont : Context LaurentNotationTest * was not created when Needs was evaluated. 
Information on the functions used can be obtained using help. 



Names [ "Laurent Functions ^ *" ] 

{AKN, AlphaKI, ANKInitialStateSetUp, BT, FiltPulse, hFiltered, InitialState, J, 
LaurentC, LaurentLK, LaurentS, M, ModulatingPulse, Modulationlndex, Modulator, 
NumberOf Curves, PhaseAngle, Receiver, ReceiverProper, S, Samplinglnterval 
StartingQuadrant, SyncSainple, T, C, ^,4/} . ' 

3 

T := 

812500 
BT := 0.3 



1 

ModulationXndex : = — 
2 

<< ModulatorData.m; 



RandomBi t S e q 

{1, 1, -1, -1, -1, 1/1, -1, 1, -1, 1, -1, 1, -1, -1, ,1, 1, 1^ 1^ _i _i ^ ^ ^ 

-1, 1, -1, -1, 1, -1, -1, 1, 1, 1, 1, -1, 1, .1, 1, 1, ,1, 1^ i/i^ { .i'/' ^' ^' 

1, 1, 1, -1, 1, -1, -1, 1, 1, 1, -1, -1, 1, 1, 1, 1, 1, 1, 1, 1, 1^ 1^ 1^ 

M-l 

SnT + AT = V £ J^^-"' CK,n'T^AT 

toiQ2 = Modulator [L] [RandomBitSeq, NumberOf Curves -►2, ModulatingPulse FiltPulse] ; 
ListPlot [{Re[tom2] , Im[tom2]} //Transpose, PlotJoined -> True] 




- Graphics - 



RandomBitSeq 

{1, 1, -1, -1, -1, 1, 1, -1/1, -1, 1, -1, 1, -1, -1, -1, 1, 1, 1, ,1, 1^ 1 1 

-1, 1, -1, -1, 1, -1, -1, 1, 1, 1, 1, -1, 1, -1, 1, _i, 1. .1, 1^ 1^ 1^ 1 .1 

1, -1, 1, -1, 1, -1, 1, 1, -1, -1, 1, -1, 1, -1, 1, 1, 1, 1, -1, -i; 1 _i, 
1, 1, 1, -1, 1, -1,. -1, 1, 1, 1, -1, -1, 1, 1, 1, 1, 1, 1, 1, 1, 1^ .1 1 _i 1 
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Receiv©r[L] ttom2, StartingQuadrant -> 4] 

{1, 1, 1, 1, 1, 1, -1, -1, -1, 1, 1, -1, 1, .1, 1^ 1^ _i _1 _i 1 1, 1, -1 -1 

1, 1, 1, -1, 1, -1, 1, -1, 1, -1, 1, 1, -1, 1, 1^ 1^ 1^ 1 ^ 

-1, -1, 1, -1, 1, 1, 1. -1, 1, -1, -1, 1, 1, 1, .1, 1^ 1^ 1^ 1^ 1^ 1^ 1^ 1^ 1^ 

We write all the possible terms 

j*o,H-5 ( Jf""-**'*"-^*"^-' *«» ) l>ulse[0] [ <5T] + 

j(a«.4*aB.j*a„.2* a«.i) i>ulsetO][T + ST] + - 4 a« ) pulse[01[2T + 6T] + 

j<OB.4*a„.3) Pulse[0][3T + <5T] + J*""-* Pulse[0][4T + <5T] + Puls©[0][5T + ST] + 
j(-«B.5) Puls©[0]t6T + <5T] + jC-as-s-a„.*) pulse[0][7T + ST] + 

j(-aN.5-aK.« -aH.7) pulse[0] [8 T + ST] + jt-H - 4 ♦ "h - 3 • OM-2*«H-i *-H Pulse[l] [ ST] + 

J(«N.4*«B-3 * «H-2 - - a ) Pulse[l] [T + <5T ] + J<-«-4--B-3 * -M - 2 * -H - 3 > Pulse[l] [2 T + ST] + 

j<«».4*«M-» -«H.4 ) Pulse[l][3T + <5T] + 5 ) Puls©[l][4T + ST] + 

j(-«i»-«> Pulse[l][5T + <5T] + jt-H-5*«H-7) Pulse[l][6T + <5T] ) 

( j(a«.**aBo*a„.a *a„ ) pulse[0][6T] + 

j(a«.4*aM.3*a«.2* p^lSefOj [.y ^ ^ J<-»- Pulse [0] [2 T + 6T] + 

j(a»-**««-3) pulse[0][3T + <5T] + J""-* Puls©[0][4T + ST] + Pulse[0][5T + ST] + 
j«-««-s) Pulse[0][6T + 5T] + jt"""-* " « ' Pulso[0][7T + <5T] + 

jC-as.s-OH-* -a„.7) Pulse[0][8T + 6T3 + j(«»-4 ♦ «»- 3 + as. a * a^. i * - ^ ) pulse [1] [ <5T] + 
j(a„.4*a„., . a^., *a„., - a^.^ ) pulse[l] [T + <5T] + 
j(a»-4*a„., * a„.3*a„.,) p^lsefl] [2 T + 6T] + 

j(aB-4*a«., -a„.,) puise[l] [3 T + ST] + jf«».*-«»-s) pulse[l] [4 T + ST] + 
j(-a„.o Pulse [1] [5 T + <5T] + * ««-7 ) Pulse[l][6 T + dT] ) / . { a„ -i^ xO, 

au.i xl, ttBj.a x2, a^-s x3, aj,.* k4, a^-s x5, an-c x6, as-? -+ x7} 

jx0.xl.x2.x3.x4 pulse[0] [<5T] + jxl-x2.x3.x4 pulse[0] [T + <5T] + 
jx2.x3.x4 Pulse [0] [2T + <5T] + J^3.x4 Pulse [0] [3T + <5T] + J^^ puise[0] [4T + (5T] + 
Pulse [0] [5T + <5T] + J-^^^ Pulse [0] [6T + <5T] + j-^5-x6 puise[0] [7 T + 5T] + 

J-X5-X6-X7 pulse[0] [8T + <5T] + jx0.x2.x3.x4 pulse[l] [6T] + jxl.x3.x4 puise[l] [T + C5T] + 
jx2.2x3.x4 puise[ll [2 T + 6T] + J^^ Pulse[l] [3 T + ^T] + J^^^-^s pulse[l] [4 T + <5T] + 
J-""^ Pulsefl] [5T + <5T] + J-^5.x7 Pulse [1] [6T + 5T] 

Clear [Absolut: e Value] 

AbsoluteValue[Pulse_] [{xO_, xl_, x2_, x3_, x4_, x5_, x6_, x7_} ] [ST 1 : = 
Module [{al, a2, a3, a4}, 
al = J«0.:^.x2.«3.x4 Pulse [0] [ST] + 

j«i.x2.x3.»4 puise[0] [T + 5T] + j-2*-3.x4 puis©[0] [2 T + 6T] + J*^*** Pulse[0] [3 T + 6T] + 
J** Pulse [0] [4 T + <5T] + Pulse [0] [5 T + ST] + J'^ Pulse [0] [6 T + *5T] + 
J-X5-X6 Pulse [0] [7 T + ST] + j-«5-x6-x7 Pulse [0] [8 T + ST] + j«o.xa.x3.x4 Pulse [1] [ST] + 
j«i.x3.xi puise[l] [T + 6T] + j«a.2x3.x4 puis©[l3 [2T + <5T] + Pulse[l] [3 T + <5T] + 
JX4-X5 Pulse [1] [4T + 6T] + J-''^ Pulse [1] [5 T + <5T] + j-'^**'' Pulse [1] [6T + <5T]; 

a2 = { Re[al] *2 + Im[al] *2 // ComplexExpand) /. {Iin[x_] :> 0, Re [x_] :> x) ; 
a2 /. x_?NuinberQy_ -> Round [x] y] 

AbsoluteValueInterfearingregi6ns[Pulse_3 [ 

{xO_, xl_, x2_, x3„, x4_, x5_, x6_, x7_> ] [<5T_] : = 
Module [{al, a2, a3, a4}, 
al = j«o.xi.x2.x3.x4 Pulse [0] [ST] + j«3*x3.x4 puise[0] [2T + 6T] Pulse [0] [6T+6T] + 

J-X5-X6-X7 Pulse [0] [8T + I5T] + j«i*-3*»4 Pulse [1] [T + <5T] + Pulse [1] [5T + <5T]; 
a2 = ( Re[al]*2 + Iin[al] ^ 2 // ComplexExpand) /. {Im[x_] :> 0, Re [x_] :> x); 
a2 / . x_ ? NumberQ y_ - > Round [x] y] 
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AbsoluteValue[Pulse[8]3 [{-1, -Ir -1, -1, -1, -1, -1, -1}][0.5T] 

x0[4]2 -HX0[12]= -2x0 [4] x0[20] +x0[20]2 -2x0 [12] x0[28] +x0[28)2 +2x0 [4] x0[36] - 
2x0[20] xO[36] +x0[36i^ +2x0[12] x0[44] -2x0[28] x0[44] +x0[44]2 -2x0[4] x0[52] + 
2x0[20] x0[52] - 2 x0[3 6] x0[52] +x0[52]2 -2x0 [12] x0[60] + 2 xO [28] xO [ 60 ] - 
2x0[44] x0[60] -hx0[60]2 +2x0[4] x0[68] -2x0[20] x0[68] +2x0[36] x0[68] - 
2 x0[52] x0[68] +x0[68]2 +2x0 [12] xl [4] - 2 xO [28] xl [4] +2x0 [44] xl[4] - 
2x0[60] xl[4] +xl[4]2 - 2 x0[4] xl[12] +2x0 [20] xl [ 12 ] - 2 xO [3 6 ] xl [ 12 ] + 
2 x0[52] xl[12] - 2 x0[e8-] xl[12] +xl[12]2 +2x0 [12] xl [20] - 2 xO [28] xl [20] + 
2 x0[44] xl[2 0] - 2 x0[60] xl[20] + 2 xl [4 ] xl [20 ] +xl[20]2 +2x0 [4] xl[28] - 
2x0 [20] xl[28] +2x0 [36] xl[28] -2x0 [52] xl[28] +2x0 [68] xl[28] -2x1 [12] xl[28] + 
xl[28]2 +2x0 [12] xl[36] -2x0 [28] xl[36] +2x0 [44] xl[36] -2x0 [60] xl[36] + 
2xl[4] xl[3 6] +2x1 [2 0] xl[36] + xl [3 6 ] ^ - 2 xO [ 4 ] xl[44) +2x0 [2 0] xir44] - 
2x0 [36] xl[44] +2x0 [52] xl[44] -2x0 [68] xl[44] +2x1 [12] xl[44] - 
2x1 [28] xl[44] +xl[44]2 + 2 xO [ 12 ] xl [ 52 ] -2x0 [28] xl[52] +2x0 [44] xl[52] - 
2 x0[60) xl[52] +2xl[4] xl[52] + 2 xl [20 ] xl [ 52 ] + 2 xl [3 6 ] xl [ 52 ] +xl[52]2 

We set T = 1 to avoid precision problems. 

T := 1 

Pulse Sampling- r= T/8 

var = Join[Table[xO[i] , {i, 0, LaurentLK[L] [0] T / Pulse Sampling -1}], 
Table [xl[i], {i, 0, LaurentLK[L] [1] T / PulseSampling - 1}]]; 
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initsol = Join[Table[I,aurentC[8] [0] [i PulseSampling] , 
{i, 0, LaurentLK[L] [0] T / PulseSampling -1}], Table [ 

LaurentC[8] [1] [i PulseSampling], {i, 0, LaurentLK[L] [1] T / PulseSampling -1}] 

$Aborted 



hessian = IdentityMatrix [Length [var] ] ; 

Pulse [8] [0] [y_] := xOLRound[y / PulseSampling] Im- 
pulse [8] [1] [y_] := xl [Round [y/ PulseSampling] ] ; 

TempFun[x_] [y_] := ( (AbsoluteValue [Pulse [8] ] [x] [y] - 1) * 2) // Expand? 

TempPun2 [x_] [y_3 : = (AbsoluteValuelnterfearingregions [Pulse [8] ] [x] [y] ) ; 

TempFun[{-l, -1, -1, -1, -1, -l, -i, -1}][0.5T]; 

ConvertToBitSeq[Deptli_] [Num_] : = 
Module [{bit, n}, 
n = Num; 

ConvertNext := (bit = Mod[n, 2]; n = Floor [n/ 2]; bit) ; 
( Table [ConvertNext, {i, 1/ Deptb) ] //Reverse) /. {0 -1}]; 

PossSeq = Table [i, {i, 1, 2^8}] // Map [ConvertToBitSeq[8] , #]&; 

We calculate the 



CostGenerator [Texi^Pun_] : = 
Module [(Cost), 
Cost = 0; 
PulseValue =0; 
CostTTpdate := (Cost - 

Cost + Map [TempPun[#] [PulseValue] &, PossSeq] // Apply [Plus, #]&//#/256 &; 
PulseValue = PulseValue + PulseSampling) ; 
Table [CostUpdate, {i, 1, T / PulseSampling // Round} ] ; Cost] ; 



AmplitudeCost = CostGenerator [TempFun] * Pulse Sait^ling; 
BERCost = CostGenerator [TempFun2] * PulseSampling; 

Save["AmplitudeCost.m", {BERCost, AmplitudeCost, PulseSampling, T}] 

In this section we find the bandwidth of the signal 
Derivative of the Sync function is represented by 

Clear [T, PulseSampling] 
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Sin[x] 
D[ , xj 

, X 

Cos [x] Sin[x] 



X x-= 
DSinc [0.0] := o 

Cos[x] Sin[x] 



DSinc [x_] /; X != 0 := 
<< AmplitudeCost .m; 

The bandwidth is given by 

TempIntO [in_, n_] : = 

L 2 



( Pulseslpling ) »^=tegrate [DSinc [ ^^^J^^^^^ (x - » P^lseSan^ling) ] 

"^""f Pulseslpling - ^l^Saaplins)], {x, 0, 9 T}] 

Teiaplntl[m_, n_] : = 

f"^r; 1. T' ] NIntegrate [DSinc [ ^ (x - m PulseSampling) 1 

V PulseSampling } ^ ^ PulseSampling wLj.»«oi«apxixig; j 

r ^ 

DSinc I ; (X - n PulseSai^pling) 1 , 

^ PulseSampling 'J 

{X, 0, 7T}] 
Table [Templnt0[0, n] , {n, 0, 5}] 

General: : stop : 

Further output of NIntegrate: :ncvb will be suppressed during this calculation. 

{1.27879x10^ 5.06325, 10.6271, 9.78983, 4.51144, -2.28746} 

PulseSampling 

1_ 
8 

ApproxBandWidth = Sum[(x0[i + 1] - xO [i] 2 / PulseSampling, {i, 0, 70}] + 
Sum[(xl[i + 1] - xlti]) *2/ PulseSampling, {i, 0, 54}]; 

Bandwidth = Sum[x0[i] x0[j] Templnt0[i, j], {i, 0, 71}, {J, 0, 71}] + 
Sum[xl[i] xl[j] Templntl[i, j], {i, 0, 55}, {j, 0, 55}]; 

Save["BandWidthCpst.m", {Bandwidth, T, } ] 
BandWidthCost = (ApproxBandWidth - 0.6) * 2; 
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BERWeight =0.3 
BandWidthWeight =0.4; 

AmplitudeWeight = 1.0 - BandWidtliWeight - BERWeight; 

0.3 

TotalCost = BandWidthWeight BandWidthCost + AmplitudeWeight AmplitudeCost + 
BERWeight * BERCost; 

dir = Directory [] ; ^ " 

SetDirectory[" . . /new2 satellite" ] ; 
<< DavidonFletcherPowell .m 
SetDirectory[dir] ; 

L := 8 ■ 

DavidonFletcherPowell [TotalCost, var, hessian, initsol, 10]; 

{0. 0207507, {DavidonFletcherPowell ^Private ^alpha$39228-> 2 . 77153}} 
{0. 0164 912, {DavidonFletcherPowell ^Private ^alpha$3 9228^ 2 .45972}} 
{0. 0142842, {DavidonFletcherPowell* Private ^alpha$39228-> 1.88521}} 
{0.0126876, {DavidonFletcherPowell* Private 'alpha$39228-» 2 ,08508}} 
{0. 0117343, {DavidonFletcherPowell * Private *alpha$39228-* 1.92086}} 
{0. 0112982, {DavidonFletcherPowell* Private *alpha$39228-» 1.55191} } 
{0. 0110916, {DavidonFletcherPowell * Private *alpha$39228-» 1. 44356} } 
{0. 010974, {DavidonFletcherPowell * Private *alpha$39228-» 1. 63222} } 
{0. 010886, {DavidonFletcherPowell * Private *alpha$39228^ 2. 1399} } 
{0 . 0108318, {DavidonFletcherPowell * Private *alpha$3 922 8-> 2 . 14114} } 
tempo [t_] = LaurentC[8] [0] [t] ; 

templ[t_] = LaurentC[8] [1] [t] ; ' 
initsolC = 

Join [Table [tempo [i PulseSampling] , {i, 0, LaurentI.K[L] [0] T / Pulse Sampling - 1} 
Table [tempi [i PulseSampling], {i, 0, LaurentIiK[L] [1] T / PulseSampling - 1} ]] ; 

ListPlot[%278[[2]], PlotJoined -> True, PlotRange -> All] 
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Graphics - 
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ListPlot[%278[[2]] //Take[#, 72]&, PlotJoined -> True, PlotRange -> All] 




0.2 - 



0.1 



10 20 30 40 50 60 70 

- Graphics - 

ListPlot [9&278C [2] ] //Take[#, -56]&, PlotJoined -> True, PlotRange -> All] 
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- Graphics - 



{AmplitudeCost, BERCost, BandWidthCost, ApproxBandWidth) /. Tliread[var -> initsol] 

{0.0687008, 0.00375044, 0.0126982, 0.712686} 

{AmplitudeCost, BERCost, BandWidthCost, ApproxBandWidth) /. Tliread[var -> Sfe278[[2]]] 

{0.0208114, 0.00872422, 0.00492789, 0.670199} 
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(ApproxBandWidtli - 0.6) ^2 /. Thread [var -> initsol] 

0.0126982 

(ApproxBandWidth - 0.6) *2 / . Tliread[var -> %278[[2]]] 

0.00492789 

OptPulse[8] [0] = Interp'6Tation[ 

Table [{PulseSamplingi, xO[i]}, {i, 0, LaurentLK[L] [0] T / PulseSampling -1)1 /. 
Thread[var -> H278[[2]]]] 

InterpolatingFunction[ { {0, 8.875}}, <>] 

OptPulse[8] [1] = Interpolation [ 

Table [{PulseSamplingi, xl[i]}, {i, 0, LaurentLK[L] [1] T / PulseSampling - 1} ] /. 
Thread [var -> ^27 8[[2]]]] 

InterpolatingFunction[ { {0 , 6.875}}, <>] 

tom3 = Modulator [L] [RandomBitSeq, NumberOf Curves -► 2, ModulatingPulse OptPulse] , 

InterpolatingFunction: idmval : Input value {-32"} ^^^s outside the range of data 
in the interpolating function. Extrapolation will be used. 

285 

InterpolatingFunction: :dmval : Input value {-32"} lies outside the range of data 
in the interpolating function. Extrapolation will be used. 

InterpolatingFunction: idmval : Input value {^^^] lies outside the range of data 

in the interpolating function. Extrapolation will be used. 
General: : stop : 

Further output of InterpolatingFunction: : dmval will be suppressed during this calculation. 
ListPlot [{Re[tom3] , Im[tom3]} //Transpose, PlotJoined -> True, AspectRatio -> 1] 




- Graphics - 

{AmplitudeCost, BERCost, BandWidtbCost, ApproxBandWidth} /. Thread [var -> initsol] 

{0.0687008, 0.00375044, 0.0126982, 0.712686} 

{AmplitudeCost, BERCost, BandWidthCost, ApproxBandWidth} /. Thread [var -> %278[[23 3 ] 

{0.0208114, 0,00872422, 0.00492789, 0.670199} 
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Plot [9&298 [t] , {t, 0, 7T}, PlotRange -> All] 

InterpolatingFunction: idmval : Input value {6.97422} lies outside the range of 
data in the interpolating function. Extrapolation will be used. 

InterpolatingFunction: :dinval : Input value {6.90039} lies outside the range of 
data in the interpolating function. Extrapolation will be used. 

InterpolatingFunction: idmval : Input value {6.93624} lies outside the range of 
data in the interpolating function. Extrapolation will be used. 

General: : stop : ^ ^* ' 

Further output of InterpolatingFunction: : dmval will be suppressed during this calculation. 



0.2 




0.05 



- Graphics - 
? Modulator 

Modulator [L] [BitSeq, Opts] assumes the following default options 

StartingQuadrant 0, InitialState ^ Table [1, {i, 1, 20} ], Samplinglnterval -> T/32, 
NumberOf Curves ^ 4 , ModulatingPulse-> LaurentC . The Pulse is assumed to have 
the following structure Pulse [L] [K] [t] 

? Interpolation 

Interpolation[data] constructs an InterpolatingFunction object which represents an 
approximate function that interpolates the data. The data can have the forms 
{{xl, fl}, {x2, f2}, ... } or {fl, f2, ... }, where in the second case, the 
xi are taken to have values 1, 2, ... . 

Table [ 

Thread [var -> %247[[2]]] 

Clear [Pulse] 

Sine [x_] /;x==0:=l 

Sinc[x_] /; X != 0 := Sin[x] / (x) 



Pulse Sampling 



T/ 8 



Pulse [8] [0] [y_] := xO[ Round [y / PulseSampling] ]; 
Pulse[8] [1] [y_] := xl [Roxmd [y / PulseSampling] ] ; 

ModulationValue[FiltPulse[8]] [{-1, -1, -1, -1, -1, -i, _i, -1}][0.5T] // N 

-0.754091 1 



ConvertToBitSeg[Depth_] [Num_] : = 
Module [{bit, n}, 
n = Num; 

ConvertNext := (bit = Mod[n, 2]/n = Floor [n/ 2] ; bit) ; 
( Table [ConvertNext, {i, 1, Depth}] //Reverse) /. {0 -> -1}]; 
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ConvertToBitSeq[8] [2] 

{-1, -1, -1, -1, -1, -1, 1, -1} 



Table [i, {i, 1, 2-^8}] // Map [Convert ToBitSeq [8 ] , #]&; 



Select [Table [i, {i, 1, 2'^8}] // Map [ Convert ToBitSeq [8] , (#[[5]] 

{{-1, -1, -1, -1, 1, -1}, 
{-1, -1, -1, -1, 1, -1, -1, 1}, {-1, -1, -1, -1/1, -1, 1, -1}, 
{-1, -1, -1, -1, 1, -1, 1, 1}, {-1, -1, -1, -1, 1, 1, -1, -1}, 
{-1, -1, -1, -1, 1, 1, -1, 1}, {-1, -1, -1, -1, 1, 1, 1, -1}, 
{-1, -1, -1, -1, 1, 1, 1, 1}, {-1, -1, -1, 1, 1, -1, -1, -1}, 
{-1, -1, -1, 1, 1, -1, -1, 1}, {-1, -1, -1, 1, 1, -1, 1, -1}, 

{-1, -1, -1, 1, 1, -1, 1, 1}, {-1, -1, -1, 1, 1, 1, -1, -1}, {-1, -1, -1, 1, 1, 1, 1}, 
{-1, -1, -1, 1, 1, 1, 1, -1}, {-1, -1, -1, 1, 1, 1, 1, 1}, 1^ X, -1, -1, -1}, 

{-1, -1, 1, -1, 1, -1, -1, 1}, {-1, -1, 1, -1, 1, -1,1, -1}, 

{-1, -1, 1, -1, 1, -1, 1, 1}, {-1, -1, 1, -1, 1, 1, -1, -1}, {-1, -1, 1, -1, 1, 1, 1}, 
{-1, -1, 1, -1, 1, 1, 1, -1}, {-1, -1, 1, -1, 1, 1, 1, 1}, {-1, -1, 1, 1, 1, _i, -1, 
{-1, -1, 1, 1, 1, -1, -1, 1}, {-1, -1, 1, 1, 1, -1, 1, -1}, {-1, -1, 1, 1, 1, -1, 1, 1}^ 
{-1, -1, 1, 1, 1, 1, -1, -1}, {-1, -1, 1, 1, 1, 1, -1, 1}, {-1, -1, 1, 1, 1, 1, 1, 
{-1, -1, 1, 1, 1, 1, 1, 1}, {-1, 1, -1, -1, 1, -1, -1, -1}, {-1, 1, -1, -1, 1, _i, -1, 1}, 
{-1, 1, -1, -1, 1, -1, 1, -1}, {-1, 1, -1, -1, 1, -1, 1, 1}, 

{-1, 1, -1, -1, 1, 1, -1, -1}, {-1, 1, -1, -1, 1, 1, -1, l}/{-l, 1, -1, -1, 1, 1, 1, 
{-1, 1, -1, -1, 1, 1, 1, 1}, {-1, 1, -1, 1, 1, -1, -1, -1}, {-1, 1, -1, 1, 1, -1, _i, 1}, 
{-1, 1, -1,1, 1, -1, 1, -1}, {-1, 1, -1, 1, 1, -1, 1, 1}, {-1, 1, -1, 1, 1, 1, -1, -1}, 
{-1, 1, -1, 1, 1, 1, -1, 1}, {-1, 1, -1, 1, 1, 1, 1, -1}, {-1, 1, -1, 1, 1, 1, 1, 1}, 
{-1, 1, 1, -1, 1, -1, -1, -1}, {-1, 1, 1, -1, 1, -1, -1, 1}, {-1, 1, 1, -1, 1, -1, 1, -1}^ 
{-1, 1, 1, -1, 1, -1, 1, 1}, {-1, 1, 1, -1, 1, 1, -1, -1}, {-1, 1, 1, -1, 1, 1, 1}, 
{-1, 1, 1, -1, 1, 1, 1, -1}, {-1, 1, 1, -1, 1, 1, 1, 1}, {-1, 1, 1, 1, 1, -1, -1, -1}, 
{-1, 1, 1, 1, 1, -1, -1, 1}, {-1, 1, 1, 1, 1, -1, 1, -1}, {-1, 1, 1, 1, 1, -1, 1, 1), 
{-1, 1, 1, 1, 1, 1, -1, -1}, {-1, 1, 1, 1, 1, 1, -1, 1}, {-1, 1, 1, 1, 1, 1, 1, -1}, 
{-1, 1, 1, 1, 1, 1, 1, 1}, {1, -1, -1, -1, 1, -1, -1, -1}, {1, -1, -1, -1, 1, -1, -1, 1}, 

{1, -1, -1, -1, 1, -1, 1, -1}, {1, -1, -1, -1, 1, -1, 1, 1}; 

{1, -1, -1, -1, 1, 1, -1, -1}, {1, -1, -1, -1, 1, 1, -1, 1}, {1, -1, -1, -1, 1, 1, 1, 
{1, -1, -1, -1, 1, 1, 1, 1}, {1, -1, -1, 1, 1, -1, -1, -1}, {1, -1, -1, 1, 1, -1, _i, 
{1, -1, -1, 1, 1, -1, 1, -1}, {1, -1, -1, 1, 1, -1, 1, 1}, {1, -1, -1, 1, 1, 1, -1}, 
{1, -1, -1, 1, 1, 1, -1, 1}, {1, -1, -1, 1, 1, 1, 1, -1}, {1, -1, -1, 1, 1, 1, 1, 1}, 
{1, -1, 1, -1, 1, -1, -1, -1}, {1, -1, 1, -1, 1, -1, -1, 1}, {1, -1, 1, -1, 1, -1, 1, _i}^ 
{1, -1, 1, -1, 1, -1, 1, 1}, {1, -1, 1, -1, 1, 1, -1, -1}, {1, -1, 1, -1, 1, 1, -1, 1}, 
{1, -1, 1, -1, 1, 1, 1, -1}, {1, -1, 1, -1, 1, 1, 1, 1}, {1, -1, 1, 1, 1, -1, -1, -1), 
{1, -1, 1, 1, 1, -1, -1, 1}, {1, -1, 1, 1, 1, -1, 1, -1}, {1, -1, 1, 1, 1, -1, 1, 1}, 
{1, -1, 1, 1, 1, 1, -1, -1}, {1, -1, 1, 1, 1, 1, -1, 1}, {1, -1, 1, 1, 1, 1, 1, 
{1, -1, 1, 1, 1, 1, 1, 1}, {1, 1, -1, -1, 1, -1, -1, -1}, {1, 1, -1, -1, 1, -1, -1, 1}, 
{1, 1, -1, -1, 1, -1, 1, -1}, {1, 1, -1, -1, 1, -1, 1, 1}, {1, 1, -1, -1, 1, 1, -1, -1), 
{1, 1, -1, -1, 1, 1, -1, 1}, {1, 1, -1, -1, 1, 1, 1, -1}, {1, 1, -1, -1, 1, a, 1, 1}, 
{1, 1, -1, 1, 1, -1, -1, -1}, {1, 1, -1, 1, i; -1, -1, 1}, {1, 1, -1, 1, 1, -1, 1, -1}, 
{1, 1, -1, 1, 1, -1, 1, 1}, {1, 1, -1, 1, 1, 1, -1, -1}, {1, 1, -1, 1, 1, 1, -1, 1}, 
{1, 1, -1, 1, 1, 1, 1, -1}, {1, 1, -1, 1, 1, 1, 1, 1}, {1, 1, 1, -1, 1, -1, -1, -1}, 
{1, 1, 1, -1, 1, -1, -1, 1}, {1, 1, 1, -1, 1, -1, 1, -1}, {1, 1, 1, -1, 1, -1, 1, 1}, 
{1, 1, 1, -1, 1, 1, -1, -1}, {1, 1, 1,. -1, 1, 1, -1, 1}, {1, 1, 1, -1, 1, 1, 1, -1}, 
{1, 1, 1, -1, 1, 1, 1, 1}, {1, 1, 1, 1, 1, -1, -1, -1}, {1, 1, 1, 1, 1, -1, -1, 1}, 
{1, 1, 1, 1, 1, -1, 1, -1}, {1, 1, 1, 1, 1, -1, 1, 1}, {1, 1, 1, 1, 1, 1, -1, -1}, 
{1, 1, 1, 1, 1, 1, -1, 1}, {1, 1, 1, 1, 1, 1, 1, -1}, {1, 1, 1, 1, 1, 1, 1/1}} 
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xl = Map[ModulationValue[FiltPulse[8] ] [#] [0.5T]&, 

Select [Table [i, {i, 1, 2*8)] // Map [Convert ToBitSeq [8] , #]&, (#[[5] ] == 1)&] ] 



{0.749266, 0.74922, 0.749173, 0.749219, 0.711482, 0.711529, 0.711482, 0.711435, 
0.829796, 0.82975, 0.829703, 0.829749, 0.792012, 0.792059, 0.792012, 0.791965 
0.776885, 0.776839, 0.776791, 0.776838, 0.739101, 0.739147, 0.7391, 0.739054, 
0.802178, 0.802131, 0.802084, 0.80213, 0.764394, 0.76444, 0.764393, 0.764347 
0.759521, 0.759475, 0.759428, 0.759474, 0.721738, 0.721784, 0.721737, 0.72169, 
0.819541, 0.819495, 0.819447, 0.819494, 0.781757, 0.781804, 0.781756, 0 78171 
0.78714, 0.787094, 0.7a7.047, 0.787093, 0.749356, 0.749403, 0.749355, 0.749309! 
0.791922, 0.791876, 0.791829, 0.791875, 0.754138, 0.754185, 0.754138, 0.754091 
0.749266, 0.74922, 0.749173, 0.749219, 0.711482, 0.711529, 0.711482, 0.711435, 
0.829796, 0.82975, 0.829703, 0.829749, 0.792012, 0.792059, 0.792012, 0.791965, 
0.776885, 0.776839, 0.776791, 0.776838, 0.739101, 0.739147, 0.7391, 0.739054, 
0.802178, 0.802131, 0.802084, 0.80213, 0.764394, 0.76444, 0.764393, 0-.764347, 
0.759521, 0.759475, 0.759428, 0:759474, 0.721738, 0.721784, 0.721737, 0.72169, 
0.819541, 0.819495, 0.819447, 0.819494, 0.781757, 0.781804, 0.781756, 0.78171, 
0.78714, 0.787094, 0.787047, 0.787093, 0.749356, 0.749403, 0.749355, 0.749309, 
0.791922, 0.791876, 0.791829, 0.791875, 0.754138, 0.754185, 0.754138, 0.754091} 

Limit [D[ Sin [x] /x, x], x -> 0] 



0 



? Limit 



Limit [expr, x->x0] finds the limiting value of expr when x approaches xO . 



Map[ 1/2 - a / 2Bt£[ — ]\&, xl] // Apply[Plus, #]& // #/128&; 




Ber[a_] := Evaluate [5&4 8] 



Ber [O.Ol] 



0.495 



Erf [Infinity] 



1 



D[Erf [x] , X] 




?Erf 



Erf [2] gives the error function erf(z). Erf[zO, zl] gives the generalized error 
function erf(zl) - erf{20). 
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Plot [Erf [X], {x, 0, 10}] 



0.9999 
0 .9998 
0.9997 
0.9996 

- Graphics - 
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?FiltPulse 

Global 'FiltPulse ■ 

T := 3 / 812500 

Null 

?7 T 

This is the symbol period 
7? T 

This is the symbol period 
T : = 3 / 812500 

Clear[T, xO, xl, x2, x3, x4, x5, x6, x7] 



'I 
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Beg^inPackage [ "riaurentrunctions * " ] 
. T::usage = "This is the symbol period" 
BT: : usage = "This is the usual product" 
h::usage = "This is the raw gaussian pulse" 
S::usage = "This denotes -modulation index Pi" 
^(r: : usage = "^/r[L,t] is Laurents function" 

hFiltered: : usage = "This is the filtered gaussian pulse" 

PhaseAngle: : usage = "This function takes some time to calculate. The following 

code will draw a graph of the functionPhaseFunction[t_] =N[ ^i. t ] ; 
phasepoints =Table [ Ct , PhaseFunction[t T] } //N, {t , 0, L ,1/40}]; 
ListPlot [phasepoints. Plot Joined -> Tarue] " 

S : :usage = "S- = Sin[ff] " 

J: : usage = " J = * " 
C::usage= "C = Cos [ff ] ■ 

M::usage="M = 2^"^" 

Modulationlndex: : usage = "ModulationJndex = h" 

Laurents :: usage = "LaurentS [L] [n] [t] = Sin[ ^[L,t + n T] ] / S" 

LaurentC: :usage = "LaurentC[L] [K] [t] is Laurents C^^t " 

AlphaKI: : usage = "AlphaKI [LL] [K, i] is Laurents Or,!" 

LaurentLK: I usage =" Laurent LK[L] [K] gives the support of CK,t " 

The start of Modulator Definitions 

AMKXnitialSt ate Setup: : usage = 
"ANKInitialStateSetUp[L] [K] [InitBitSeq, AccumulatedPhase] sets up the sequence 
of prior states of Ak,m that the modulator went thgough to get to the 
constelation point specified by AccumulatedPhase which is really Ao,o" 

AKN: :usage = "AKN[L] [K] [{State, AccumulatedPhase}] defines Ak,^ in terms of Ao,h" 

ModulatingPulse: :usage = 
"The Pulse is assumed to have the following structure Pulse [L] [K] [t] " 

NumberOfCurves: : usage = "The number of pulses used by the modulator" 

Samplinglnterval: : usage = 
"Samplinglnterval is the interval between samples of the output of the modulator" 

InitialState: rusage = "InitialState is the set of bits that are assumed 
to be present before i.e {a_i,a.2# -••}" 

StartingQuadrant: : usage - "StartingQuadrant = Ao,.iand is a number" 

Modulator: :usage = "Modulator[L] [BitSeq,Opts] assumes the following 
default options StartingQuadrant 0, InitialState Table[l, {i, 
1,20}] , Samplinglnterval -> T/32, JSTumberOf Curves 4 , ModulatingPulse-*^ 
LaurentC . The Pulse is assumed to have the following structure Pulse [L] [K] [t] " 

Start of the Receiver Functions 
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PiltPulse: : usage = -The default pulse and is called by PiltPulse [L] [K] [t]- 

SyncSample: :usage = "Given tliat the sampling interval is T/32,then sync 

sample has rang -16 to 16 and this moves the point used to demodulate " 

Receiver: : usage = ** Receiver [L] [Input Seq, Opts] assumes the following" 
default options { Start ingQuadrant-^0, InitialState-^ {1, 1, 1, 1, 1, 1, 1, l, 
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1} , Samplinglnterval -> T/32,ModulatingPulse 
FiltPulse, SyncSample -> O^NumberOf Curves -> 2 " 

ReceiverProper: lusage - "ReceiverProper[L] [ StartingQuadrant, Initiaistate, 
San^lingXnterval^ModulatingPulse, SyncSample,NumberOf Curves, Input Seq] 
is called by Receiver after all the options have been resolved** 

Begin [" ^ Private ^ ** ] 

VLog[2] 

a : = 

2 TT BT 



h[t_] : = 



V2 TT a T 



Ideally we would and did define the effect of the convolution of h[t] by the formula below. However, this version of 
Mathematica gives an error. 

T • 

lifil [t_] := Release [Module [{r}, J"^ ^^^^ ^^]] 

hFiltered[PulseWidth_] [t_] := Module [ {xl, x2, x3}, xl = 

/Krr*i ^m^xr-ni^i T/..-. f^l^Ctl-r] . PulseWidth PulseWidth T 
(N[#l, 40] &) [Table [{tl, J ^ dr], {tl, , _ , "20"'^^^' 

x2 = Interpolation [xl] ; x2 [t] ] 

S :=N[ModulationIndexjr] 

C := Cos [»] ; 
S := Sin[«] ; 
J (e^» //Chop); 
M := 2^-^; 

PhaseAngle [L_] [t_] /; t ^ 0 : = 0 
PhaseAngle [L_] [t_] /; t ^ I*T := » 

PhaseAngle[I._] [t_] : = 

PhaseAngle [I.] [t_] =Module[{xl, x2, x3, x4, x5, x6>, xl = hPiltered[3 I, T] [ tl - -^1; 

2 

Evaluate [xl] dtl}, {t2, 0, L T, ]1; Interpolation [x2] [t]! 

^^^tl*-* t_] /;0<t<LT:= PhaseAngle [L] [t] 

t_] /; 2 LT > t ^ LT := » - PhaseAngle [L] [t - LT] 

We need to put this to avoid the function being extrapolated 

^^r[L_ , t_] /; ! (0 < t < LT) && ! (2 I.T > t i LT) : = 0 
Laurent S[L_] [n_] [t_] := Sin[ ^[I., t + n T] ] / S 
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AlphaKI [LL_] [K_, i_] /; (0 < i < LL) && (0 <= K < 2^'^) : = 
Module[{xl, x2, x3, KNum) , xl : = ^x2 = Mod[KNu2n, 2] ; KNum = ™^ " x2j ; 

KNum = K; x3 = Table [xl, {ii, 0, LL - 1} ] ; x3 ] 

LaurentLK[L_3 [K_] : = 
Module [{xl}, xl = Table [L {2 - AlphaKI [L] [K, ii] ) - ix, {ii, 1, L-l)]; Min[xl]] 

LaurentC[L_3 [K_] [t.] /J^ 0' <= K < 2^ : = 

L-l 

LaurentS[L] [0] [t] LaurentS[L] [ii + L AlphaKI [L] [K, ii] ] [t] 

The start of the modulator function 

ANKInitialStateSetUp[L_] [K_] [InitBitSe<i_, Acc\imulatedPhase_] : = 
Module [{xl, x2, x3, x4, x5, acuphase, initbitseq} , 
initbitseq = InitBitSeq; 
acuphase = AccvunulatedPhase; 
UpdateSeq : = 

Module[{}, xl = acuphase - Sum [initbitseq [ [i ] ] AlphaKI [L] [K, i] , {i, 1, l -1}]. 
acuphase = acuphase - First [initbitseq] ; initbitseq = Rest [ initbitseq]; xl]'- 
Table [ XJpdateSeq, {i, 1, LaurentLK[L] [K] } ] ] 

AKN[L_] [K_] [{State_, AccuinulatedPhase_} ] : = 
AccximulatedPhase - Sum [State [ [i + 1] ] AlphaKI [L] [K, i] , {i, 1, L -1}] 

Options [Modulator] := {StartingQuadrant 0, InitialState Table [1, {i, 1, 20}], 
Samplinglnterval T/32, NumberOf Curves 4, ModulatingPulse LaurentC} 

Modulator [L_] [Bit Seq_, Opts ] : = 

Module [ 

{xl, x2, x3, x4, x5, x6, state, AccumalatedPhase, seq, AKNState, Curves, Pulse} 
xl = Samplinglnterval/* {Opts} /. Options [Modulator] ; ' 
state = InitialState /• {Opts} /. Options [Modulator] ; 
x3 = StartingQuadrant/. {Opts} /. Options [Modulator] ; 
x4 = Samplinglnterval/. {Opts} /. Options [Modulator] ; 
Pulse = ModulatingPulse /. {Opts} /• Options [Modulator] ; 
Curves = (NumberOf Curves / . {Opts} /. Options [Modulator] ) - 1; 
AccumulatedPhase = x3 ; 
seq = BitSeq; 
Table [ 

AKNState[K] = ANKInitialStateSetUp[L] [K] [state, AccumulatedPhase], {K, 0, Curves}] - 
x5 := Module [{}, state = Join [ {First [seq] } , Drop [state, -1]]; 
AccumulatedPhase = AccumulatedPhase + First [seq]; 
seq = Rest [seq]; 

Table [AKNState [K] = Join [ {AKN [L] [K] [ {state, AccumulatedPhase} ]} , 

Drop [AKNState [K] , -1]], {K, 0, Curves}]; 
x6[r„] = Sum[Sum[ ( j) A«»staf[K] ni * H ] pulse[L] [K] [r + i T] , 

{i, 0, LaurentLK[L] [K] -1}], {K, 0, Curves}]; 
Table [x6[r], {r, 0, T - x4, x4}]]; 
Table[x5, {kk , 1, Length [BitSeq] } ] //Flatten] 

Options [Receiver] := {StartingQuadrant 0, 

InitialState-* {1, 1, 1, 1, 1, 1, 1, i, i, i, i, i, i, i, i, i, Samplinglnter^ 
ModulatmgPulse -+ FiltPulse, SyncSample -> 0, NumberOfCurves -> 2}; 

Receiver [L_] [InputSeq_, Opts ] : = 

Module [{xl, x2, x3, x4, x5, x6}, 

xl = StartingQuadrant/. {Opts} /. Options [Receiver] ; 
x2 = InitialState/. {Opts} /. Options [Receiver] ; 
x3 = Samplinglnterval/. {Opts} /. Options [Receiver] ; 
x4 = ModulatingPulse / . {Opts} /. Options [Receiver] ; 
x5 = SyncSample/- {Opts} /. Options [Receiver] ; 
x6 = NumberOf Curves / . {Opts} /. Options [Receiver] ; 
ReceiverProper[L] [xl, x2, x3, x4, x5, x6. Input Seq] ] 
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ReceiverProper[L_] [ Start ingQuadrant_, InitialStato_, Samplinglnterval_, 
HodulatingPulse_, SyncSainple_, NumberOf Curve s_, IiiputSeq_] : = 

Module [{xl, sgn, ReceivedSeq, ExpectedValue, seq, ReceiveNext, 3D, ^} , 
xl = T / Samplinglnterval; 

ReceivedSeq = Partition [Input Seq, xl] // Transpose //#[ [SyncSample +l]]t; 
ExpectedValue = 

ModulatingPulsetL] [0] [ (Laurent LK[L] [0] /2) T + SyncSample Samplinglnterval] 
Si = StartingQuadrant; 
sgn =0; 

^ = O; 

seq = ReceivedSeq; 
ReceiveNext := 
Module [{xl, x2}, 

xl = ( (-1)"*^ First [seq]) // Im; 
If[Abs[xl - ExpectedValue ] <= Abs [xl + ExpectedValue ] , 

ID = Join[ B, {1}]; ^ = + 1, D = Join[ E), {-!}]; ^ = ^ -1 ]; 
seq = Rest [seq]; sgn = Mod[sgn + 1, 2]]; 
Table [ReceiveNext, {i, 1, Length [ReceivedSeq] }] ; 
D] 

End[] 

EndPackage [] 



